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SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 


This report summarizes the results of one phase of a study of techniques for 
measuring atmospheric ozone that was carried out at NASA Wallops Flight Center. This 
study represents the second phase of a program designed to improve techniques for the 
measurement of atmospheric ozone. This phase of the program was concerned with two 
distinct tasks: 

1) A study of the sensitivity of Dobson direct sun measurements and the ozone 
amounts inferred from those measurements, to variations in the atmospheric 
temperature profile. The study used the plane-parallel, Monte-Carlo model 
developed and tested under the initial phase of this program (1) and a 
series of standard atmospheres taken from Reference 2. 

2) Development and testing of a Monte-Carlo simulation of radiative transfer 
through a spherical shell atmospheric model. Spherical atmospheric models 
are necessary to study the Dobson Umkehr measurements which use large solar 
zenith angles. 

Under task 1 a plane-parallel Monte-Carlo model of the Dobson measurements systems 
developed under the initial phase of this program was used to study the effect of 
variations of atmospheric temperature profiles on the total ozone inferred from Dobson 



direct sun measurements. The ozone profile and total ozone in each atmospheric model 
were held constant. The results of the study are summarized in Table 1. The errors 
listed are for vanishingly small instrumental field of view, which essentially minimizes 
errors due to single and multiple scattering. The errors listed then reflect the 
sensitivity of the ozone amounts, inferred from direct sun measurements, to errors in the 
assumed mean ozone absorption coefficients, and, to a lesser extent, errors in the 
molecular and aerosol extinction. 


TABLE 1 

ERRORS IN OZONE AMOUNT INFERRED USING DOBSON GUAGE 
(OVERESTIMATE = POSITIVE) 


Lines 

Std. 

Atm. 

Winter 

Warm 

Subarctic 
Winter 
Cold Jan 

July 

Midlatitudes 
Jan July 

Subtropical 
Jan July 

Tropical 



60°N 

60°N 

60°N 

60°N 

45°N 

45°N 

30°N 

30°N 

1 5°N 

A 

-1.49 

-2.23 

-3.90 

-3.19 

-1.54 

-2.74 

-1.74 

-2.36 

-2.05 

-2.25 

C 

-0.33 

-0.72 

-3.77 

-2.46 

0.24 

-1.69 

-0.14 

-1.18 

-0.74 

-1.10 

D 

0.75 

-1.43 

-2.55 

-2.09 

-0.41 

-1.23 

-0.23 

-0.32 

-0.27 

-0.45 

A-D 

-2.13 

-2.56 

-4.30 

-3.59 

-1.93 

-3.18 

-2.23 

-2.96 

-2.63 

-2.84 

C-D 

-0.01 

-0.19 

-4.71 

-2.73 

0.75 

-2.06 

-0.10 

-1 .84 

-1.11 

-1.61 

A-D 

-3.05 

-3.49 

-3.82 

-3.69 

-3.04 

-3.54 

-3.04 

-3.30 

-3.10 

-3.16 


A source of the errors in the ozone estimations is the temperature dependence of 
the ozone absorption coefficient. Essentially the ozone absorption coefficient used in 
reducing the direct sun data is an average over the ozone density distribution in the 
atmosphere. Since the ozone absorption coefficient varies with temperature and hence 
with height in the atmosphere we will arrive at a different average value for different 
temperature profiles or ozone profiles. 

Under task 2 a Monte-Carlo simulation of radiative transfer through a spherical ly- 
symmetric, layered atmosphere was developed and tested against well known solutions. 

That part of the study will be described in a separate report and is not covered here. 

Recommendations: It is recommended that future work should emphasize the following 

tasks: 

1) The ozone absorption spectra of ozone in the Hartley band system (2900A to 

o 

3400A) should be investigated as a function of wavelength and temperature to 
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determine whether or not a wide band line pair would reduce or eliminate the 
ozone absorption temperature dependence. 

2) As an interim approach the appropriate values of the average ozone coefficient 
a should be evaluated for every standard temperature model for a number of 
possible ozone profiles; from these investigations appropriate values of A^ 
and to be used in Dobson reductions should be derived. 


INTRODUCTION 


The amount of ozone present in the earth's atmosphere is of considerable interest to 
the scientific community. Among the reasons for the importance of ozone is its role in 
the transfer of solar energy in the atmosphere, absorbing harmful ultra-violet radiation, 
and providing a source of thermal energy in the stratosphere. Ozone is also important 
as a tracer used to study atmospheric circulation. 

Ozone measurements have been made for a number of years using a variety of instru- 
ments and techniques. The Dobson spectrophotometer has served as a standard ground- 
based instrument for the measurement of total atmospheric ozone for over two decades. 

The Dobson gauge is a double-monochromator that isolates two narrow bands in the ultra- 
violet spectrum. The wavelengths of the bands are chosen so that one wavelength is 
strongly absorbed by ozone while the other, longer wavelength is affected very little 
by ozone absorption. Measurements of the intensities reaching the instrument are used 
to calculate the amount of ozone that the solar radiation has passed through in arriving 
at the instrument. 

In order to assess the accuracy of the total ozone inferred from Dobson measurements 
a study was initiated at Wallops Flight Center several years ago. The first phase of 
that study was concerned with the effects of single and multiple scattering by molecules 
and particles in the atmosphere on the signals received by the Dobson gauge. That 
phase was completed last year and concluded that single and multiple scattering can 
contribute significant errors to the amount of ozone inferred from direct sun measurements. 
However, these errors could be reduced to a minimum if not completely eliminated by using 
a narrow instrument field-of-view and using line pair coupling. We also found a signifi- 
cant sensitivity, of the ozone amount inferred, to atmospheric temperature profiles. 

This report describes the second phase of the study which was concerned with the errors 
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in the total ozone amount estimated from Dobson direct sun measurements for standard 
atmospheres ranging from tropical through sub-arctic conditions. 


DOBSON DIRECT SUN MEASUREMENTS 


System Description 

The Dobson spectrophotometer is a specialized double beam monochromator used for 
determining atmospheric ozone by measuring the ratio of the intensities of ultra-violet 
light at two wavelengths in the solar spectrum. The wavelength pair is selected so that 
one wavelength is much more strongly absorbed by ozone than the other, so that the 
intensity ratio can be used to estimate the total amount of ozone in the optical path 
from the sun to the instrument. By appropriate adjustments of the spectrophotometer 
one of several specified wavelength pairs can be selected for the measurement, depending 
on illumination conditions, to provide an estimate of the amount of atmospheric ozone. 

The Dobson instrument can be used in three distinct modes: a direct sun measurement, 
a "zenith sky" measurement and an "Umkehr" mode. In the direct sun mode a prism or "Sun 
Director" is used to project the sun's image onto the entrance window of the Dobson 
spectrophotometer. The direct sun mode is the preferred method for measuring total ozorie. 
Details of the Dobson instrument and its use in all three modes are given in Reference 1 
and will not be repeated here. In this report we are only concerned with direct sun 
measurements in a plane parallel atmosphere. 

The Dobson instrument measures the relative intensities of the two wavelengths in a 
given pair by first separating them in a spectroscope, isolating them with slits, and 
attenuating the beam at the wavelength least absorbed by ozone with a variable thickness 
optical filter (the filter wedge) . The two beams are brought back into coincidence with 
a second spectroscope so that they impinge on a photo-electric null-meter. The null- 
meter does not measure absolute intensities, instead it indicates when the difference 
between the intensities is zero. The null-meter consists of a chopper wheel in the 
optical path, which passes the two beams sequentially and a photo-multiplier upon which 
one beam and then the other fall sequentially as allowed by the chopper wheel. An 
amplifier is used for the A.C. output of the photo-multiplier, then a rectifier converts 
the alternating current to direct current, and finally a micro-ammeter displays the 
direct current. 


4 



The ratio of the intensities of the two beams is "measured" by adjusting the filter 
wedge until the micro-ammeter reads zero. Then the intensity ratio is computed from the 
position of the filter wedge and the instrument calibration data. 

Wavelength Doublet and Wavelength-Doublet Pairs . - The Dobson machine can be 
adjusted to make measurements at any of five wavelength doublets. Table 2 taken from 
Reference 3 gives the pertinent data on the wavelength doublets A, B, C, D and C\ At 
present wavelength doublet B is not used, so that only four wavelength doublets are 

used. Also, when making observations for total ozone, it is desirable to make nearly 
simultaneous observation on two wavelength doublets. The wavelength doublet pairs 
used are AD, CD, and CC'. 


TABLE 2 


TABLE OF WAVELENGTHS AND OTHER CONSTANTS 
(absorption coefficients are based on log 10 ) 


Designation 

of 

Wavelength 

doublet 

Mean 

Wave- 

length 

A.U. 

Ozone 

Absorpti on 
Coeffi cient 
(atm cm) - ' 

Atmospheric 
Scattering 
Coeffi cient 
(atm) -1 


a 

a 1 

1 

a-a 

e 

6' 

b~b' 

e-e' 

a-a 1 

Short 

3055 

1.882 


0.491 



A 



1.762 


0.116 

0*066 

Long 

3254 

0.120 


0.375 



Short 

3088 

1.287 


0.470 

* 


B 



1.223 


0.113 

0.092 

Long 

3291 

0.064 


0.357 



Short 

3114.5 

0.912 


0.453 



c 



0.865 


o.no 

0.127 

Long 

3324 

0.047 


0.343 



Short 

31 76 

0.391 


0.416 



D 



0.374 


0.104 

0.278 

Long 

3398 

0.017 


0.312 



Short 

3324 

0.047 

■H 

0.343 



C 







Long 

4536 

Nil 

■ 





5 





Theory of Measurement Using Direct Sunlight . - For a plane-parallel stratified 
atmosphere (Fig. 1) containing Rayleigh, scatte re rs (molecules), ozone molecules and 
aerosols, the change in sunlight from one layer to the next can be expressed as: 

<11 (h) = -1(h) {a ni (h)p m (h) + a o (h)p Q (h) (1) 

+ o (h)p (h)} sec 2 d£. 

r r 

The quantity 1(h) denotes the intensity at altitude h while cr(h) denotes a cross 
section for scattering or absorption in cm 2 /molecule and p(h) denotes the number of 
molecules or particles per cubic centimeter. Subscripts denote molecules (m) , ozone 
molecules (o) , and particles (p). Note that in this analysis we have neglected any 
light that might be scattered into the path from other parts of the atmosphere; we 
only considered losses from the beam due to scattering and absorption. 


LOCAL VERTICAL SUN 
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( 2 ) 

(3) 

(4) 

(5) 

In I(h s ) = In I (o ) - {e m T m + §pTp + a X} sec Z (6) 

Essentially the quantity T ffl represents the "scale height" of the atmosphere, that is, 
the thickness of the molecular atmosphere if it were reduced to a uniform layer at 
Standard Temperature and Pressure. Similarly, then represents the Rayleigh volume 
scattering coefficient for the atmosphere at Standard Temperature and Pressure. The 
volume coefficient can be calculated from the Rayleigh-Cabannes formulae for the 
scattering cross section per molecule and Loschmidt's number n Q , the number of molecules 
per cm 3 at Standard Temperature and Pressure. The product e m T m then is the normal or 
vertical optical thickness of the atmosphere due to Rayleigh or molecular scattering. 
Similarly 6 p T p represents the normal optical thickness of the atmosphere aerosol and 
a X represents the normal optical thickness of the ozone layer. The parameter X is . 
the equivalent thickness of ozone or the ozone scale height, and is the quantity of 
interest. The quantity a is the average value of the ozone absorption coefficient 
over the ozone profile. In this plane parallel model, restricted to zenith angles 
less than 60°, we do not concern ourselves with the effect of refraction on the path 
length, and we neglect the effect of the earth's curvature. 

In Dobson's derivation he sets 6 = e m T ffl and <5 = <s p T p . Then we have: 

In I(h s ) = In I { o) - (6 + 6 + a X} sec Z (7J 


Integrating the equation from the top of the sensible atmosphere (h=o) to the 

h„ 


surface (h=h g ) we have 


In 


I(h $ ) 


1 ( 0 ) 




{a m (h)p m (h) + a 0 (h)p Q (h) + a p (h)p p (h ) > sec Z dh 


If we now define the dimensionless quantities: 

h 


/” 
o /*h 

]'• 


a m (h)p m (h) dh = 6 T 
m f m m m 


P (h S (h) * • Vp 


and 


f 


s 0 (h)p 0 (h) dh = « X 


then the equation for the intensity reaching the ground based instrument is: , 
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where & now denotes the vertical optical thickness of the atmosphere due to molecular 
scattering and 6 represents the vertical optical thickness due to scattering 
by particulate matter. 

In Dobson's derivation for the direct sun measurement the quantity sec Z multiplying 
the expression for the ozone optical thickness and molecular optical thickness is 
replaced by quantities y and m. The parameter m is defined to be the equivalent path 
length through the molecular atmosphere allowing for atmospheric refraction and the 
curvature of the earth. When the zenith angle is zero, m is defined to be 1.0. The 
quantity y is the relative path length of sunlight through the ozone layer which is 
presumed to be at 22 km above the earth surface. When the zenith angle is zero, y is 
also defined to be 1.0. However in the final equations used m and y are taken to be 
essentially the same as sec Z. Since we confine ourselves to zenith angles less than 
60°, these distinctions are unnecessary. Then in the interest of brevity and clarity we 
will use only sec Z. 

The preceding derivation used natural logarithms while the original derivation of 
Dobson employed logarithms to the base ten. Changing to log 10 we have 

log 10 I x = logio I QX - (e + § + a X) sec Z. 


In this form of the equations, the coefficients are related to the preceding coefficients 
through relations of the form: 

« 10 = = 0-4342 v 6io = 0.4342 e e and 6 10 = 0.4342 6 g . 

Writing the equation for two wavelengths x and x' and subtracting, we arrive at: 

logio log! - |(e - s') + {<s - <s') + (a - a') xj sec z (9 
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We have: -N 


x 

or X 


(a - a') x X + (e - e') x + (6 - 6') x 


sec Z 


(a-ot' )xsec Z 


(g-e')x+ (6-6 ')x 


(a-a')x 


OD 

( 12 ) 


The subscript x denotes the wavelength dependence of the terms while primes indicate 
the longer wavelength of the line pair. 


The term containing (6 - <$') cannot be evaluated explicitly because the attenuation 
coefficients for particulate scattering usually are not known. As a means of avoiding 
this difficulty, observations are often made nearly simultaneously on two wavelength 
pairs, say A and D, and the resulting equations (3-11) subtracted and solved for X. 

The result for the AD line pairs, using the coefficients given in Table 2 is: 


X AD = 0.7205 (N a - N d )/ sec Z -0.009 (13) 

- |(<S-i5 i ) a - ( <5-5 ‘ ) qJ / j(ot-a 1 ) A - (a-a')^J 

The last term in equation (13) has been set to be zero, in Dobson's derivation, 
on the assumption that particulate scattering coefficients are only slightly wavelength 
dependent, thus cancelling the two terms in the numerator. We then have 


X AD = 0.7205 (N a - N d ) / sec Z - 0.009 (14) 

which is identical to Dobson's equation 5.1 in Reference 3 and is the basis for the 
reduction of AD line pair measurements. A similar equation can be derived for the line 
pairs CD and AC. 


MONTE-CARLO SIMULATION 


A Monte-Carlo model of the Dobson ozone estimation scheme was assembled and tested 
in the initial phase of this study. Details of the Monte-Carlo model are given in 
Reference 1 and will not be reported here. Briefly, the model can compute the effect of 
multiple scattered sunlight on the atmospheric ozone density inferred from Dobson spectro- 
photometer measurements. The program possesses a useful feature. Unlike the earlier 
Monte-Carlo simulations of Kattawar and Plass 4 and Collins and Wells 5 the present method 
does not use direct sampling of photons. Rather, it samples hypothetical photons from 
each scattering on a probability basis. This technique has two important advantages: 
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1) The variance is considerably reduced. 

2) We can sample sky radiance in exactly specified directions and do not have to 
use averaging techniques. 


Polarization is fully treated by applying the Mueller algebra to the Stokes vector 
associated with each photon. The program can model ten receiver configurations simulta- 
neously, Including direct sun observations from ground level, zenith sky observations 
from ground level and nadir observations from satellite altitudes for any solar zenith 
angle less than 60°. The constraint on the solar zenith angle is imposed by the plane- 
parallel atmosphere assumed in the Monte-Carlo simulation. Other receivers at any 
altitude and with any orientation can also be specified by the user. Thus, as a byproduct 
of the Dobson computations we can investigate the effect of multiple scattering on obser- 
vations made in or above the atmosphere. 


The Atmospheric Models 

The computer program used in this study will accept as data any atmospheric model 
a user wishes to Investigate. In order to limit the amount of computer time used we 
restricted ourselves to a series of Standard Atmospheres given in Reference 2. The 
models used are listed in Table 3 and the temperature and molecular density profiles 
for each model are shown in Fig. 2 through 6. The variation of the molecular 
densities of the models from the standard atmosphere density is shown in Fig. 7. 


TABLE 3 


MODEL 

ATMOSPHERES 


Designation 

Time of Year 


Standard Atmosphere 



Tropical 15°N 



Subtropical 30°N 

January 

July 

Midlatitude 45°N 

January 

July 

Subarctic 60°N 

January 

July 

Subarctic 60°N 

Winter Warm 

Winter Cold 
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ALTITUDE (KM) 








ALTITUDE (KM) 



Figure 4 - Vertical Temperature Profile for Midlatitude (45*N) 
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ALTITUDE {KM) 











Our standard atmosphere used the July temperature model of the U.S. Standard Atmosphere 
while the molecular densities were taken from Champion's atmospheric model (6) which 
represents mean conditions throughout the year for latitudes near 30 degrees above 30 km. 
Below 30 km we used an average of the Cole and Kantors (7) 30 degree latitude summer and 
winter atmospheres to provide a smooth fit between Champion's densities and those of 
the U.S. Standard Atmosphere (2), which were used exclusively below 10 kms. For the 
other models, molecular densities were calculated from tabulated values of atmospheric 
density and mean molecular weights as functions of altitude (2). For aerosol densities 
we employed the Elterman model (8). The size distribution used was assumed to be that 
of Junge; i.e., the number of particles with diameters in the range D to D + dD was: 


n(D)dD = (v-l)D 0 v_1 D v dD = D > D q 

=0 D < D„ 

0 

In this investigation we used: 


D q = 0.2u 
v = 4.0 


The aerosol density was kept constant in all of the atmospheric models investigated. 
The aerosol scattering properties are computed with a Mie scattering computer program 
supplied by Holland (9). We assumed that all aerosol particles were spherical with a 
real refractive index of (1.55-0i). Wavelength dependence of the aerosol refractive 

index was ignored so that the sole effect of changing wavelength was to alter the 

TT D 

Hie size parameter X =— for the particles in the distribution. The ozone density 
profile used in these studies was taken from data by Dutsch (10). The ozone absorption 
cross sections used were taken from reference 8 and are summarized in Table 4. The 
variation of the ozone coefficients with temperature is shown in Fig. 8. For tempera- 
tures intermediate to the published values, linear interpolation was used. The ozone 
density was kept constant for all of the atmospheric models studied. 

In these simulations we assumed that all other molecular components behaved as 
Rayleigh scatterers, with scattering cross sections per molecule as given in Table 4. 
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TABLE 4 


ABSORPTION CROSS-SECTIONS FOR OZONE 
AND RAYLEIGH SCATTERING CROSS-SECTION 
(Values based on log e ) 


21 9°K 
18°C 

LINE 

°ABS 

DOBSON PAIRS 

2 

X(y) cm /molecule 

.3055 1 .757xl0~ 19 

"A" Lines 

.3254 1.526x1 0 -20 

-?o 

.31145 8.634x10 

"C" Lines 

.3324 4. 578x1 O' 21 

.3176 3.745xlO" 20 

"D" Lines 

.3398 1 .451x1 0 -21 


229°K 
-44° C 

21 4° K 
-59°C 


0 ABS 

0 ABS 

°R 

2 

cm /molecule 

2 

cm /molecule 

2 

cm /molecule 

1.608x1 O' 1 9 

1.582xl0~ 19 

5.187xlO" 26 

1. 254x1 0" 20 

1.221xl0‘ 20 

3.964x1 0" 26 

7.778xl0" 20 

7.554xl0“ 20 

4. 798x1 0" 26 

3.461xl0‘ 21 

3.312xl0" 21 

3.621xl0" 26 

3.259xlO“ Z0 

3.223xl0 -20 

4. 396x1 O" 26 

1 . 116xl0 -21 

1 .042x1 0" 21 

3. 297x1 O' 26 


The Rayleigh scattering cross sections were calculated from the equation: 


3.327x10 

,4.258 


-28 


7 

cm /molecule where X, the wavelength is expressed in microns. 
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MOLECULE 


"A" LINES 


LEGEND 

Q A - 305SA' 

□ A * 32S4A*’ 

^ A ■= 3114. SA* 

"C" LINES 

^ A = 3324A” 

© A * 3176A" 

"D" LINES 

Q A « 3398A" 



TEMPERATURE DEGREES K 


Figure 8'- Ozone Absorption Cross Sections VS Temperature 
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Analysis of Dobson Direct Sun Measurements 

For each model atmosphere described above we simulated radiative transfer through 
the model for each wavelength in the line pairs A, C and D. The parameter defined 
in equation 9 was calculated as a function of receiver field of view using the sum 
of the direct and scattered intensities reaching the receiver in the direct sun mode. 

The parameters N x so calculated were then used to determine the total ozone in the 
atmosphere through Dobson's equation using the Dobson coefficients given in Table 2. 

That is, we calculated the quantities eqs. (13) or (14). 

x x = (sec Z)-' - [(B-8'l + («-6')] U-a 1 )" 1 

■ \ "x z >'' - B x 

Now in equation (10), N x is defined as: 

- logio 

and appears to depend on the values of the solar flux outside the sensible atmosphere for 
the three line pairs used. Normally these ratios are inferred by making observations at 
a large number of zenith angles or air masses and extrapolating to zero air mass. However 
It is easy to show that the scattered intensity for any receiver field of view is always 
directly proportional to the intensity incident at the top of the atmosphere as is the 
direct attenuated intensity so we can write I = K I Q and I 1 = K ' 1 1 Q . Hence N x = 

Togio (I 0 /I 0 ') - 1og 10 (I/I 1 ) = log 10 (K'/K). Or equivalently we can set I Q = I 0 ' for 
each wavelength pair so that = log^o (K'/K). The values of K and K 1 are calculated by 
the Monte-Carlo simulation. Thus the values N x are independent of the ratios of the 
incident intensities. 

In our calculations we know the true value of the equivalent ozone thickness X 
since we could integrate the model ozone profile over height to obtain the integrated 
column count. Dividing by Loschmidt's number, the number of molecules in a cubic 
centimeter at Standard Temperature and Pressure, gave us the value of X. We could then 
obtain the true value of a, the average ozone absorption coefficient, by performing the 
integral in eq 5 namely: 

/* h s 

a X = J 0 Q (h) P Q (h) dh 
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and divide by X to obtain a in (atm cm)”^ or cm -1 SIP. In our model the value of X 
was 0.36567 cm STP. 

In the same way we could obtain the proper values of 6 and <S, and so calculate the 
correct values of A x and B x eq. (4) used in Dobson's equation, for each model atmosphere. 
Details of the calculation of the coefficients are contained in Appendix A. 


DISCUSSION 


The results of our model simulations are summarized in Table 5. These results are 
for extremely small instrument fields of view so that errors due to singly or multiply 
scattered light are minimized. These results then reflect the sensitivity of the ozone 
amounts inferred from Dobson direct sun measurements , to atmospheric temperature profiles. 


The source of the errors can be appreciated by neglecting the effects of scattering 
by molecules and aerosols. We have shown in Reference 1 that the errors due to scattering 
can be essentially eliminated by narrowing the instrument field of view and by using line 
pair coupling. Then, viewed as a simple absorption experiment, the Dobson gauge measures 

I_ = io -olX or logio aX » for a zenit ^ sur1, O 5 ) 


Now the quantities aX are defined by equation (5) as: 

h 


»X = J" 3 Q (h) p 0 (h)dh 


(16) 


where o (h) is the ozone absorption coefficient (cr,,* 'molecule) and p (h) is the ozone 
nunber density /Mlecula) 


\ 


cm J 




Now 


/ s 

n Q X = J p Q (h)dh 


(17) 


is the integrated column count or the nunber of ozone molecules in a column with a base 
area of one square centimeter reaching to the top of the atmosphere. The symbol n 0 
denotes Loschimdt's number (2.687xl0 19 molecules/ cm 3 ) the number of molecules per cm 3 
at standard temperature and pressure. So we can write 


A 

J °o (h) p o (h)dh ~ Vo X 


«X 


(18) 
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TABLE 5 


PERCENTAGE ERROR (OVERESTIMATE) USING DOBSON COEFFICIENTS 


MODEL 

A LINES 

C LINES 

D LINES 

A-D LINES 

C-D LINES 

A-C LINES 

Standard (30°N) 

-1.49 

0.33 

0.75 

-2.13 

- 0.01 

-3.05 

45°N Jan 

-2.74 

-1 .69 

-1.23 

-3.18 

-2.06 

-3.54 

45°N July 

-1.74 

-0.14 

-0.23 

-2.23 

- 0.10 

-3.04 

60°N Winter Warm 

-2.23 

-0.72 

-1.43 

-2.56 

- 2.01 

-3.49 

60°N Winter Cold 

-3.90 

-3.77 

-2.55 

-4.30 

-4.71 

-3.82 

60°N Jan 

-3.19 

-2.46 

-2.09 

-3.59 

-2.73 

-3.69 

60°N July 

-1.54 

0.24 

-0.41 

-1 .93 

0.75 

-3.04 

30°N Jan 

-2.36 

-1.18 

-0.32 

-2.95 

-1.84 

-3.30 

30°N July 

-2.05 

-0.74 

-0.27 

-2.63 

- 1.11 

-3.10 

15°N 

-2.25 

- 1.10 

-0.45 

-2.84 

-1.61 

-3.16 

STATISTICS 







Maximum 

-1.49 

0.33 

0.75 

-1.93 

0.75 

-3.04 

Minimum 

-3.90 

-3.77 

-2.55 

-4.30 

-4.71 

-3.82 

Mean 

- 2.20 

- 1.12 

-0.82 

-2.83 

-1.54 

-3.32 

Standard Deviation 

0.76 

1.26 

1.02 

0.72 

1.58 

0.29 


Overall Mean = -2.00 


Overall Standard Deviation = 0,99 



where a Q denotes an average ozone absorption cross section averaged over the ozone 
density profile. 

But what we really have in hand is the quantity aX and to arrive at a correct value 
for X we need the correct value of a or But since a Q varies with temperature and 
so with height in the atmosphere, then in order to calculate the true average value, we 
need to know the ozone density profile with height. 


In practice the value of X we get will depend on the value 
as the value of the absorption coefficient for ozone depends on 
height. If the temperature dependence could be eliminated then 
absorption coefficient from the integral eq. (16) and have 


A / 

J 0 P 0 (h)dh = a ^P 0 


Q (h)dh = oti 0 X = aX 


of a that we use, as long 
temperature and hence on 
we could remove the 


According to Mitra (11) the temperature sensitivity of the ozone absorption 
coefficients has been Investigated by Wulf and Melvin (12) from -78°C to 250°C. They 

o o 

found that in the Hartley Bands between 2900A to 3400A when the temperature was decreased 
there was an Increase in the maxima in the bands and a decrease in the minima. Similar 
investigations by Vassy (13) for temperatures between -80° to 20°C did not show a 
corresponding increase in the band maxima, but confirmed the decrease in the minima. 

Vassy also investigated the effects of varying pressures (between 19 nrn Hg to 760 mm Hg) 
and found practically no pressure effect. The absorption data of Vigroux quoted in 
Reference 8 show a consistent decrease in both maxima and minima with decreasing tempera- 
ture. These conflicting results should be investigated in the laboratory, since Wulf 
and Melvin's results, if correct, suggest that a broad band line pair might be chosen so 
that temperature effects are eliminated or considerably reduced. However the results 
of Vassy and Vigroux appear to contradict Wulf's results. 


Table 6 compares the correct values of the parameters A^ and calculated for 
each temperature model with those of Dobson (3). The discrepancy in the values 
of A, reflect the temperature sensitivity of A = ■ # ' V/rv or , ■ J _ t 

say to the average value of the ozone absorption coefficient. Another approach to 
this problem is to study the variation in values A^ for each standard atmospheric 
temperature model for various reasonable ozone density profiles with height. In a 
sense one could derive an average value for the standard models that would yield better 
estimates of total ozone than the values now in use. 
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TABLE 6 

CORRECTED VALUES OF A^ AND 
COMPARED WITH DOBSON VALUES 


Lines 


A 

C 

D 


A-D 


C-D 


A-C 


Subarctic 



Std. 

Atm. 

Winter 

Warm 

Winter 

Cold 

Jan 

July 



60°N 

60°N 

60°N 

60°N 

A x 

0.5754 

0.5781 

0.5879 

0.5836 

0.5739 

B x 

0.0665 

0.0656 

0.0668 

0.0663 

0.0650 

A x 

1.1488 

1.1560 

1.1920 

1.1763 

1 .1448 

B x 

0.1257 

0.1242 

1.1396 

0.1263 

0.1229 

A x 

2.6885 

2.7129 

2.7441 

2.7316 

2.6843 

B x 

0.2842 

0.2822 

0.2854 

0.2841 

0.2791 

\ 

0.7322 

0.7346 

0.7482 

0.7422 

0.7300 

B x 

0.0098 

0.0096 

0.0096 

0.0096 

0.0096 

A x 

2.0016 

2.0143 

2.1074 

2.0658 

1.9961 

B x 

0.0035 

0.0034 

0.0034 

0.0034 

0.0034 

A x 

1.1513 

1.1561 

1.1601 

1.1586 

1.1509 


0.0062 

0.0061 

0.0061 

0.0061 

0.0061 



Midlatitudes Subtropical Tropical 


Jan 

July 

Jan 

July 


45°N 

45°N 

30°N 

30°N 

1 5°N 

0.5817 

0.5754 

0.5800 

0.5776 

0.5787 

0.0665 

0.0655 

0.0666 

0.0659 

0.0661 

1 .1638 

1 .1507 

1.1660 

1.1584 

1.1625 

0.1265 

0.1240 

0.1267 

0.1251 

0.1256 

2.7213 

2.6874 

2.7070 

2.6932 

2.6972 

0,2849 

0.2804 

0.2848 

0.2817 

0.2820 

0.7398 

0.7322 

0.7382 

0.7352 

0.7368 

0.0096 

0.0096 

0.0097 

0.0096 

0.0096 

2.0520 

2.0124 

2.0480 

2.0326 

2.0432 

0.0035 

0.0034 

0.0035 

0.0035 

0.0035 

1.1569 

1.1511 

1.1541 

1.1518 

1.1523 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 


Dobson 

Coefficients 
Ref (3) 

.5675 
.0660 
1 .1560 
0.1270 
2.6730 
0.2780 
0.7205 
0.0090 
, 2.0370 
0.01 20 
1.1147 
0.0064 



In summary then we have reached the following conclusions. 


1) The equations used in the reduction of Dobson direct sun measurement depend 
on the atmospheric temperature profile and the ozone density profile. 

This dependence is due to the temperature dependence of the ozone absorption 
coefficient. 

o o 

2) The absorption spectra of ozone in the Hartley band system (2900A to 3400A) 
should be investigated as a function of wavelength and temperature to 
elucidate the discrepancy in the published data, determine whether or not 

a wide band line pair would reduce or eliminate the ozone absorption 
temperature dependence and to provide more reliable data for interpolating 
between temperatures. 

3) As an interim approach, the appropriate values of the average ozone coefficient 
a should be evaluated for every standard temperature model for a number of 
possible ozone profiles; from these investigations aporopriate values of A^ 
and B a to be used in Dobson reductions should be derived. 
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APPENDIX A 


DERIVATION OF THE COEFFICIENTS OF THE TOTAL OZONE ESTIMATION EQUATION 


Referring to equation (12) we see that the form of the ozone estimation equation is 


where 


X * A. N, cos Z - B. 

A A A 

\ = logio^-i^- logio^fr^ 


(A-l) 

(A-2) 


and A, and B are constants which may be derived from the properties of the atmosphere 
for each line pair or set of coupled line pairs. Values of A^ and have been derived 
for standard conditions making certain assumptions in regard to scattering and absorption 
cross sections of the atmospheric constituents. In this section we explain how we used 
the information generated by the Monte-Carlo program "INPUT" subroutine to derive new 
values of these coefficients which would be relevant for our particular model atmospheres 
and ozone temperature sensitivity coefficients. 


Five sets of coefficients were derived for the following cases: 

1) A lines 

2) C lines 

3) D lines 

4) A-D lines 

5) C-D lines 


We shall treat the line pairs first and then explain the procedure for the coupled 
line pairs. 


Equation (12) shows that for uncoupled line pairs we may write 


and 


1 

a-a ' 


(0-0') + (0-fi') 



(A-3) 

(A-4) 


Referring back to equation (8) we can see that aX can be related to the contribution of 
absorption to the optical thickness of the atmosphere. Neglecting aerosol and molecular 
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(A-5) 


scattering we would obtain for an overhead sun (i.e., zenith angle equal to zero) 

iogi</-r~) * - « x 


But the contribution of absorption by ozone to the total optical thickness can be 
represented as 


T oz = l0 ^e’ 


Cf) ■ - 2 - 3 ’ 3 

= 2.303 «X 


(A-6) 


Hence, we can write a-a' as 


T _ _ - t 

f _ oz 

a-01 " 2.303X 


oz 


(A— 7> 


Now in the computer program we compute specifically the total optical thickness in 
each case and the contribution of scattering to this value, r $ . Since absorption in 
the wavelengths of interest is attributed to ozone alone we have that t Q2 = - t s . 


Thus ( A- 7 ) may be rewritten as 


V T s 


- w 


2.303X 


(A-8) 


All of the parameters of the right hand side of equation (A-8) are available from the 
atmospheric model and cross section assumptions and are explicitly evaluated and output 
on a routine basis. Thus we have 


(a-a') ^ = 


2.303X 

T t' T s • T t’ +T s' 


(A-9) 


By a similar argument we are able to derive the value of from the scattering contribu 
tion to the total optical thickness. If we neglect absorption by ozone in equation (8) 
we obtain for an overhead sun 



- (e+s) 


(A-10) 
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(A-ll) 


The contribution of scattering to the total optical thickness is in this case. 


Thus 


Hence, we can write 


t_ = log -2.303 log 10 y*- 

s e 1 X 4 ox 

= 2.303 (0+s) 

T s 

e+5 " '27353 

(0+6) - ( 8 1 +6 1 ) as (t s -t s '}/2.303 and (A-4 becomes) 

(t s -T ' ) 

R = ^ S - A 

C X 2.303 X 


(A-12) 


(A-13) 


where A^ is computed by equation (A-9). 


Equation (A-9) and (A-13) are employed to evaluate corrected coefficients A x and 

B for the A, C and D lines. For line pair coupling we rewrote (A-l) as 

X 

X B x 

N a cos Z - + jf- (A-l 4) 

X X 

When we couple line pairs we subtract another equation of this type from (A-16) obtaining 


W x -V> COS 



(A-l 5) 


where the dash now stands for another line pair. Inverting (A-l 5) we obtain 

x ■ ■ *7) > cos z ■ (k ~ V) (k ■ v) (A -’ 6 ’ 

Equation (A-16) gives us the appropriate coefficients for coupled line pairs in terms of 
parameters evaluated in the computer model set-up. These are employed to obtain the 
coefficients attributed to Thomas in the computer output and are applied to evaluate 
errors as a function of the field of view. 
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